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Abstract 

Neutrino and anti-neutrino states coming from the neutral current 
or Zq decay are blind with respect to the flavor. The neutrino oscil¬ 
lation is observed and formulated when its flavor is known. However, 
it has been shown that we can see neutrino oscillation pattern for 
Zq decay neutrinos provided that both neutrino and anti-neutrino are 
detected. In this paper, we restudy this oscillation via quantum field 
theory approach. Through this approach, we find that the oscillation 
pattern ceases if the distance between the detectors is larger than the 
coherence length, while both neutrino and antineutrino states may be 
coherent. Also the uncertainty of source (region of Zq decay) does not 
have any role in the coherency of neutrino and antineutrino. 


1 Introduction 

While neutrino oscillation is a window to the new physics, it is one of the 
most interesting quantum mechanical phenomena. Historically, neutrino os- 
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cillation has been established more than 50 years and confirmed experimen¬ 
tally more than 10 years [TJ. The observation of neutrino oscillation depends 
on the coherency of neutrinos during the production, propagation and detec¬ 
tion M- The production and detection coherence conditions are satisfied 
provided that the intrinsic quantum mechanical energy uncertainties during 
these processes are large compared to the energy difference A Ej^ of different 
neutrino mass eigenstates: 

Ami , . 

A- Ejk ~ 2 E ^ aE ' W 

where a# = mm{a^ od , This condition implies that, during the pro¬ 

duction and detection processes, one cannot discriminate the neutrino mass 
eigenstates. Conservation of the coherency during the propagation means 
that the wave packets describing the mass eigenstates overlap from the pro¬ 
duction until the detection regions. The wave packets describing the dif¬ 
ferent neutrino mass eigenstates propagate with different group velocities. 

After propagating L, the separation of different mass wave packets is 2E 2 L. 
Consequently, the coherent propagation is guaranteed provided that 


A m% 

2 E 2 


L O'xu 







( 2 ) 


where v g is the average group velocity of the wave packets of different neutrino 
mass eigenstates and a xu is their common effective spatial width. In other 
words, similar to the double-slit experiment, if one could determine which 
mass eigenstate is created or detected, the neutrino oscillation pattern would 
disappear. For instance, conservation of energy and momentum implies that 
exact determination of energy-momentum of charged leptons leads to the 
determination of mass eigenstate of the corresponding neutrino (in fact exact 
momentum conservation causes neutrino state is entangled kinematically to 
the corresponding charged lepton state) and neutrino oscillation is ceased [4j. 

The kinematics analysis shows that a neutrino state created through 
charged current interactions has a specific flavor. For instance, muon neutrino 
is created by the pion decay while muon decay gives only electron neutrino. 
In contrast, neutral current or Z° decay is blind with respect to the neu¬ 
trino flavors. In other words, every flavor eigenstate as well as every mass 
eigenstate is created with the same probability. However, there is another 
property that is noticeable; neutrino and antineutrino states are entirely cor¬ 
related in the sense that they have same flavor. It has been shown that if 
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both neutrino and antineutrino are detected, one can observe neutrino oscil¬ 
lation pattern between the detectors p|. Nevertheless, if only either neutrino 
or antineutrino is detected, the neutrino oscillation is ceased; therefore, it 
is a realization of the Einstein-Podolsky-Rosen paradox [6] . In the center of 
mass frame, in particular, the oscillation pattern occurs at distance L + L, 
where L and L are the distance of the neutrino and antineutrino detectors 
from the source, respectively. In this paper, we reanalyze it via quantum 
held theory approaches. In this approach, the oscillating states become in¬ 
termediate states, not directly observed, which propagate between a source 
and a detector. The localization conditions are respected with attributing 
a localized wave function to interacting initial and final states in the source 
and detector dusmuDun]. Indeed, these localizations are essential for the 
observation of the neutrino oscillation and guarantee the coherence issues 
|4j. In the case of Zq decay neutrinos, we will see the localization of source 
(region of Zq decay) is not important while the localization of detectors plays 
role in the coherency condition of neutrino and antineutrino. Moreover, as 
was said in general, the coherency is spoiled during propagation because the 
group velocities of various mass eigenstates are different. Therefore, maybe 
one expects that when both neutrino and antineutrino propagate coherently, 
it is possible to have oscillation pattern. However, we will show that it is 
necessary the distance between detectors to be smaller than the coherence 
length. 

In the following, we develop Z 0 decay neutrino oscillation through quan¬ 
tum held theory approach. Finally, we discuss on the coherency properties 
which appears through quantum held theory approach. 

2 Developing neutral current neutrino oscil¬ 
lation through quantum field theory ap¬ 
proach 

We can describe any particle physics processes by S-matrix formalism in 
quantum held theory provided that it is adjusted according to the physical 
situations. In particular, to describe neutrino oscillation one needs to notice 
that neutrinos are produced and detected in confined space-time regions. 
The source and detector regions are separated by a finite distance which is 
usually much larger than the size of these regions. 
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Neutral current neutrino oscillation consists of the following three pro¬ 
cesses: 

• Creation neutrino and antineutrino in the source 

Z 0 —* v + P. 

• Detection of neutrino in the corresponding detector 

v -f- Dj —^ D p -\-1 . 

• Detection of antineutrino in the other detector 

v Dj —y Dp + Z~C 

In order to define the initial and final states, the localizations of interactions 
in the source and detectors require to integrate on momentum with a lo¬ 
calized distribution function around the corresponding averaged momentum. 
Therefore, the initial states are defined as follows: 

I Z 0 ) = J [dp]F 2 (p, P)|Z 0 , p) , 

\Di)= j[dk]F Dl {k,K)\D I ,k) , 

\D:)= I[dk]F Dl (k,K)\Dpk) , ( 3 ) 

where D] (Dj) is target in the detectors of neutrino (antineutrino), [dp] 
denotes — d p - The final states are written similarly as follows: 

(2ir) 3 ^/2Ep J 


Dp) = f [dk']Fp) F (k', K')\D F , k') , 

Dp) = I [dk']Fp) F (k', K.')\D f , k') , 

1 1~) = j[dk' l ]F l -(k",K")\rX) , 

|Z + ) = I [dk”]Fi+ (k", K") |Z + , k") . (4) 
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Here, D F (Dp) refers to created nucleon (antinucleon) in detector due to 
neutrino (antineutrino) collision. l~ (Z + ) denotes the created charged lepton 
corresponding to neutrino (antineutrinos) in the detector. In the above de¬ 
fined states, F’s are momentum distribution functions which are localized 
around the corresponding mean momentum. The amplitude of the neutrino 
and anti-neutrino production - propagation - detection processes is given by 
the following matrix element: 

iA a g = (D F ,l-,D F ,l + \fe~ l f d4 ^)_ i\ Di ,D I: Z 0 ) 

= {D F ,r,D F ,l + \i J d A x x j d A x J d A x 2 Ui{x 2 )Ui{x)Ui{x 1 )\D I) D F Z 0 ) 

= ^ £ KPuPf, r, Dp, D It Zo) . (5) 


where T is the time ordering operator and HLi are the weak interaction Hamil¬ 
tonian. The quantity A p ' w ' is the plane wave amplitude of the process with 
the j’th neutrino and antineutrino mass eigenstate propagating between the 
source and the detectors and is written as follows: 


iA p j - w ik,k',k",k,k', k") 


x 

x 

X 


i / d A x x / d A x / d A x 2 M jD (k,k ', jfc")e- i (*-*'-*"> sa 


d A q fi + 


rrin 


(27r) 4 q 2 — m 2 + ie 


e -iq {x2 - x ) M ^ z ^ e - ipx 


f A i + m J 

J (27r) 4 q 2 — m 2 + ie 


Mj D (k, k', k")e 


—i(k—k'—k")x i 
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( 6 ) 


where M’s are the plane wave amplitudes of the processes. It is convenient 
to switch to shifted 4-coordinate variables x, x x and x 2 defined according to 


X\ —» X\ + Xp> X 2 ^X 2 + X F) X ^ X + Xp, 

where the propagation times T and T are dehned by 


tp ~ tp — T, 
and the propagation distances by 
x fl -xp = L, 
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xp - x P = L, 




and we redefine 


F[k) = f(k)e ikXD ==► f(k) = F{k)e~ ikXD . 


Taking into account that 

fi + rrij = }-m j = J^Vj(q, s)vj(q, s), 


we redefine the amplitudes (including spinors) correspond to the production 
and neutrino and antineutrino detection processes, respectively, as follows: 


M jjZ (p) 


Vj(q) 

\/2<fo 


M z (p) 


Uj(q) 

VWo' 


M jD (k, k', k") = l -^±M D (k, k ', A:"), 

v 2go 

M jD {k,k'X) = M D {k,k'X) 1 ^. 

V Zc lo 

Substituting the initial and final state from (J3J) and (0} into (JSJ) and using 
above issues, we have 


A 


a/3 


X 

X 



A 2 g0 -UtoT+i g.L 

(27t) 4 g 2 — m 2 + ie 


f A _ 2 ^0 -tq 0 r+tq.L 

J (27t) 4 g 2 — m 2 + ie 





d A x\e~ iqxi 


d A x 2 e ~ iqx2 


x J [dk\f Dl (k,K)e- ikxi J [dk'U^XX')^ J [dk!']XXX"y k " xi 

x J [dk]f Dl (k, K)e~ ikX2 j [dk']r DF (k', K r )e ik ' X2 J [dk"]/; + (k", K")e^ 
x M jD {kX,k")M jjZ {p)M jD {k,k',k") . (7) 


Notice that the integration over x in the recent equation leads to the 5- 
Dirac function representing energy-momentum conservation in the source. 
Hereafter, we assume, for simplicity, the momentum wave functions of the 
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initial and final states to be Gaussian which are sharply peaked around the 
corresponding averaged momentum similar to 


/(P,P*) = 




3/2 


-(p-pp 


(Try 


( 8 ) 


where cr p , width of momentum distribution, is assumed to be very smaller 
than the corresponding averaged momentum. Therefore, similar to the method 
presented in [7j, the amplitude of the total process can be written as 


AaP ^ U a j Upj J 


d A q 


2go 


(27r) 4 q 2 — m 2 + ie 


3 — iqoT+iq.L, 


X 


X 


d A q 


2Qo 


D -iqoT+ici.L 


(27t) 4 q 2 — m 2 + ie 
d 4 x i exp [-f (g 0 + E Di 


E 


o. 


[dp]fz(p, P)(27t) 5 (q + q-p ) 
^ z -)t 1 + i(q + K-K , -K")x 1 


(xi - v^fi ) 2 (xi - v Di? fi ) 2 (xi - v,-fi ) 5 








x 


I d 4 x 2 exp [—f(g 0 + E Di - - E l+ )t 2 + i(q + K — K' — K")x 2 

(x 2 - v Dl t 2 ) 2 (x 2 - vp F f 2 ) 2 (x 2 - v z +f 2 ) 2 




4cr 2 n 

xDp 


Ki + 


x M jS (fc, F, k")M jjZ (p)M jD (k, F, k" 


( 9 ) 


where a x ’s are the position uncertainties which are related to the momentum 
ones through a x cr p ~ - and v’s denote the group velocities of the correspond¬ 
ing particles. Since the elements of matrix M are smooth functions of the 
on-shell 4-momenta, whereas the wave packets of the external states are as¬ 
sumed to be sharply peaked at or near the corresponding mean momentum, 
one can replace M by their values at the mean momenta and pull out of the 
integral. Moreover, we define 


Ed! ~ E Df — Ei - E 


Di 


E d — E d — Ei+ — En 


Di 


K - K' - K" = K 


Di 


K —K' — K" = K S , 
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.2 / V Dj V Df Vi- 

_2 ' _2 ' 2 /’ 


_ 2 / V -Dj I 

= ^d(-2— + “2— T ~2 

°xD, °xD f °x) 


_ 2 ( ' Ui 

V P = <7 xD (- 2 


v 5/ v 5 f Vh_x 
' _2 ' _2 /’ 


V >/ | V Pf | V Z 

' _2 ' 


— a xP 1^2 1 „2 _2 

^ — 


V?_ 


V '2 _z _z 

TcP, °xP F °xZ 


® xDi ® xDp ^xl+ 

V _ _2 / V L V 5 f V j+ x 

S D = + ^2~ + ^2“)- 

^sDj a xP F xZ+ 


- ± p - 

Therefore, using above issues and carrying out the integration over x\ and 
X 2 one can write the amplitude as follows: 

1 


A a g OC 


X 


V3 


5^ U* aj u 0 j M jD (K, K\ K")M jjZ (P)M jD (K, K', K")(2n) 
j 


d i( l _ 2 ^° -iqoT+iq.L _ 

,/ (27r) 4 q 2 — m 2 + ie 


(27r) 4 q 2 — m 2 + ie 


d A( i _ 2( ?0 c -iq 0 T+iq.L 


x /'[dp]/.(p,P)h 4 (g + g-p)e 5(g) e s(,?) 


( 10 ) 


where 


and 


S( (J ) = ( K ^ + q ) 2 + K®> + ^)- I K D + q | v D ] 2 


4rr 2 

^ a pD 


^Id X d 


_ (Kp + q ) 2 ^ [(go + E d )— | K d + q | Vp ] 2 


Kd 


4cr 2 f=,An 


pD 


D 


( 11 ) 


( 12 ) 


with \fj(D) = S D ( 5 ) —v 2 D 0 y Now, one should carry out the integration over 
the momentum of either propagating neutrino or propagating antineutrino, 
ffere, we integrate over the momentum of antineutrino. After applying the 
following change in integration variable 


p - q = p 


p = p 


d 6 p = d 3 p\ 


we have 


A a/3 cx -= U^M^K, K\ K")M jjZ (P)M jD (K, K\ K")( 2 'k) a 


x 


j 

Pq 


2 go 


(27t) 4 g 2 — m 2 + ie 


_ e -i qoT+ i^ e -s{ q ) Lp']/ Z ( p ' + q p) e SV) 


x- 


2 Ep> 


j—iEjT+ip'.l. 


p ' 2 — rrZ + ie 


(13) 






















where Ej = E p — q 0l p ' 2 = E 2 — p'\We perform the integral over p' using the 
Grimus-Stockinger theorem |5J 


d 3 p' 


0(p')e 


zp'.L 


, L —^oo 


27r 2 


P'a ~ P ' 2 + ie 


0(Pa L ) 




where L = This theorem is valid for a function ej) which is differentiable 
at least three times such that (j) itself and its first and second derivatives 

decrease at least as \ as 

p u 


p' |—> oo. Performing the recent stage, one can 


write the amplitude as follows: 
2vr 2 1 


A 


o3 


OC 


X 


E KjU Pj M jD (K, K 1 , K")M jjZ (P)M jD (K , K'. K")( 2vr) 4 
v j 

d \ 4 f {Ep ~ go) e -Wor+iq.L e -S( g ) /z( p . + q)e -5fe) e -*|pdL e -^f 


xe 


(27t) 4 g 2 — m 2 + ze 

iq 0 T 


(14) 


in which | pj |= y E 2 — rrij 2 . To carry out the integration over the neutrino 

4-momentum, it will be more convenient for us to integrate first over q° 
and then over the components of q. It is noticeable that since the pole at 
q° = —Ej + ie is not physical, the contribution to the integral is only given 
by the residue at the pole of the neutrino propagator at q° = Ej — ie . We 
obtain 


A a f) oc 


4f7T 2 1 


E KjUpiMjDiK, K\ K" ) Mjj Z {P)Mj d (K, A", K")(2n) 


x e lEpT 


d 3 q 
(27r) s 


fz(pj + q )(E p - Ej{ q))e-^( q )T-r) el q . L e -i|p,|Z 


xe~ s ^e~ s< ^ Pj \ 


(15) 


The remaining integration over q can be done by using saddle-point approx¬ 
imation at q = Pj. We expand Ej( q) about q = p ; as follows: 

E M) = Ej(Pj) + (q - P j)vj + ( 16 ) 

where Vj = at | q |=| p ? |. Also S'(q) + S(p — q) is expanded as follows; 
5(q) + 5(p - q) = S(pj) + S(pj) + ^ 9 ^ S) ( q _ Pj .)2 + ( 17 ) 
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where the first derivative of S+S at q = Pj vanishes and the second derivative 
is given by 


d 2 (S + S) _ 1 (vj-vp) 2 1 ( v j — v _d) 2 _ n 

9q 2 2 a 2 pD 2a 2 D \ D + 2a 2 pD + 2a 2 pD X D r { ’ 

Using the above issues, one can perform the integration over d 3 q. Conse¬ 
quently, the amplitude is obtained as follows: 

Aity 2 1 - - - 

A,s oc —-=^(/;,(/ Sj Af jC (A',A",A'")M„ z (P)M, D (A',A'',A'")(2 I r)‘‘ 


xe iEpT 


x exp 


—iPL 


(E p -Ej)f z ( P)e 


-iEj [T — X 1 ) + fpj(U + E) 


((L+Z)- Vj -(T-T )) 5 


21L 


x exp (—S(pj) - Stpj))- 


(19) 


The probability of the process is proportional to | ^4. Q( g 2 . In a practical ex¬ 
perimental setting L and L are usually fixed and known quantity while T and 
T are not measured. Therefore, the probability of detecting a neutrino with 
flavor a and an antineutrino with flavor (3 by the neutrino and antineutrino 
detectors located at the distances L and L from the source, respectively, is 
obtained by the time average of | A a p | 2 , which leads to 


p°e « *)(i + £)-(sfe) + sfe)) 


j,k 


■( Sip k ) + S(p k )) 


i(L + L)(Ej - E k )(Q k Vj + OjVfc) 


(flfcv 2 + QjV 2 ) 

(L + T) 2 (vj - v fc ) 2 (Ej - E k ) 2 ^ k ■ 
2(O fc v 2 + Mjvl) 2(fl fc v 2 + 


2 \ I ’ 


( 20 ) 


where 


4 ?7 r 2 _ _ - — 

Nj = -j-M jB (K,K')M jjZ (P)M jD (K,K')(27T)E- iEpT (E p -E j )e- i ^ L f z (P) 


and N? is its complex conjugate. Since we are concerned with relativistic 
neutrinos, we use the following approximations. The differences between 
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the energies and momenta of various mass eigenstates are due to the thin 
splitting of masses. Hence, we approximate 


Ej — E 




P 


2 E 1 


( 21 ) 


in which E is the common neutrino energy when rrii = 0 and p is determined 
from the energy-momentum conservation [T2]. Equation (12T|) leads to the 
following approximations 




Pj-E + (p-l)^, 


and 


Vj - 1 




2 E 2 ' 


Also, due to these approximations, one can easily show that 


( 22 ) 

(23) 


a i - 


(24) 


where a 2 = a 2 D + a 2 q and 


uj = 1 + 


v ^) 5 




+ 


a 


xD 


(i 


vp ) 2 


D 




D 


(25) 


Moreover, one can see that the relativistic approximation leads S + S to 
be minimum. Therefore, in the relativistic approximation we obtain the 
following expression for the flavor-changing probability: 


Pa/3 


5>< 

i,k 


a] 


UfjU ak Ul k NjNl exp 



(26) 


with the oscillation length L° s k c 
given by 


T OSC 

L jk 


AirE 

H’ 


and the coherence length L ™ k , 


r coh _ 

^jk — 


,— ‘IE 2 

2 V ^\A^-\ a - 


for j ^ k, 


(27) 


The exponent in the transition probability obtained for neutral current neu¬ 
trino includes three terms; the first term leads to the usual oscillation pattern 
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between the detectors, the second term indicates that the coherency condi¬ 
tion is satisfied provided that the distance between the detectors is not larger 
than the coherence length and finally the third term shows that the position 
uncertainty due to the detection mechanisms must not be larger than the 
oscillation length. It is noticeable that 

• the coherent propagation of both neutrino and antineutrino is not suf¬ 
ficient because the oscillation pattern is ceased if the distance between 
the detectors is larger than the coherence length. In fact, in quantum 
field theory approach, the conservation of energy-momentum due to the 
integration over the coordinates of the Z 0 decay vertex makes neutrino 
and antineutrino propagators entirely entangled. 

• the integration over the coordinates of the vertex of Z 0 decay gives 
energy-momentum conversation and the uncertainty of source is, prac¬ 
tically, excluded from calculations. In other words, the source uncer¬ 
tainty does not play any role in the coherency of neutral current neutri¬ 
nos and the detector uncertainties are analogues to the production and 
detection uncertainty in the case of he standard neutrino oscillation in 
the baseline L + L. 


3 Summary and Discussion 

It has been shown that we can see neutrino oscillation pattern for Z 0 decay 
neutrinos provided that both neutrino and anti-neutrino are detected [5]. In 
this paper, we restudy this oscillation and corresponding decoherence issues 
via quantum field theory approach. We should emphasis that although, de¬ 
tection of two neutrinos is far from the experiment, the theoretical study of 
the neutral current neutrino oscillation leads to some nontrivial viewpoints 
about the theory of neutrino oscillation. In quantum field theory approach, 
neutrino and antineutrino are described by free propagators and the initial 
and final particle states are described by corresponding wave functions. The 
conservation of energy-momentum due to the integration over the coordi¬ 
nates of vertex of the Z 0 decay makes neutrino and antineutrino propagators 
entirely entangled. Therefore, the coherency of individual neutrino and an¬ 
tineutrino is not enough for oscillation, but the distance between the corre¬ 
sponding detectors have to be smaller than the coherence length. The other 
important result is related to the uncertainties of source and detectors; the 
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source uncertainty does not play any role in the coherency of neutral current 
neutrinos and the detector uncertainties are analogues to the production and 
detection uncertainty in the case of he standard neutrino oscillation in the 
baseline L + L. 
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